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Abstract 

We find that with uniform mesh, the numerical schemes derived from finite 
element method can keep a preserved symplectic structure in one-dimensional case 
and a preserved multisymplectic structure in two-dimentional case in certain discrete 
version respectively. These results are in fact the intrinsic reason that the numerical 
experiments indicate that such finite element algorithms are accurate in practice. 
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As both the finite elemente method and symplectic scheme || as well as 
multisymplectic scheme are powerful tools to solve differential equations numericaly, 
it is interesting to explore if there is some relation between them. We will give a partial 
answer to this question by a simple example. 

In order to show the symplectic or multisymplectic structures in the scheme derived 
from finite elemente method , we consider the boundary value problem of the semi- 
linear elliptic equation in one- dimensional and two-dimensional spaces: 

Ati = f(u) in Q, u\dn = on dfl. (1) 

where Q is a bounded domain in ]R n , n = 1,2 and f(u) is nonlinear and sufficiently 
smooth enough function. The weak formulation of the boundary value problem of the 
equation is to find u : Q —>■ Hq(Q) such that 

[ Vu-Vvdx = - [ f -vdx Vt>G#o(ft), (2) 
Jn Jn 

Let 

a(w, v) = I Vu-Vvdx, (f,v) = — f ■ vdx, 
Jn Jn 

then (0) becomes, 

a{u,v) = (./». (3) 

It is important to note that the equation is in fact an ODE with Lagrangian on IR 1 or 
Lagrangian PDE in IR 2 respectively. The details of the symplectic and multisymplectic 
structures in the Lagrangian formalism can be found in @ || @] || || ||. 

In one-dimensional case we first discrete M 1 with regular lattice E} with equal 
spatial step h 

]L ={•••, Xi-i, Xi, Xi+i, ■■ ■} 

Let O be a segment in IR 1 and (fi linear shape function such that (fi(xj) = 6ij. As 
usual, Ui = u(xi). 

From the finite element method, we get the scheme 

Uj+i - 2m + Uj-i f + 1 £/ 1 ^ 
h 

The right side can be rewritten as 



/ f(y\uktpk)tpidx. (4) 
Ji - 1 i-i 



h ■= / f(ui-i(fi-i + Uiipi)(pidx + / f(ui<pi + u i+1 (pi+i)(pidx. (5) 

Ji—l Ji 

It should be noted that the variables u^s in the scheme can be released from the 
solution space to the function space by means of relevant discrete Euler-Lagrange 
(DEL) cohomologically equivalent relation |4j]. Therefore, as long as working with the 
DEL cohomology class associated with the DEL euation (times by certain 1-form) the 
Mfc's can be regarded as in the function space in general rather than in the solution 
space. 

Introducing the DEL 1-forms 

E Di := {u i+ i - 2ui + Ui-i - hlijdui, (6) 



such that the null DEL 1-form gives rise to the equation in the finite element method. 
The DEL condition reads 



dE Di = 0. (7) 

Namely, the DEL 1-forms are closed. It is straightforward to see that from the DEL 
condition it follows 

dv,i + i A dui + du^i A du^ 
= h([ /'(f/i-i^i-i + u i Lpi)Lpi- 1 Lpidx)du i -i A dui 

Ji-l 
ri+1 

+ h( f'{uifi + Ui + i<f i+ i)<fi + i(pidx)du i+ i A dui, (8) 

Ji 



i.e. 



where 



u^=J l \ (9) 



ri+l 

) _ (i _ j f'(ui<pi + u i+1 <pi + i)<pi + i(pidx)diH+i A du^ (10) 

It can be checked that this 2-form is closed w.r.t. d on the function space and its 
coefficients are non-degenerate, so that it is a symplectic structure for the scheme 
derived from finite element method in one-dimensional and it is preserved. 
In the two-dimensional case, the semi-linear equation becomes 

u XlX1 + u X2X2 = f(u) in Q, u\ 9n = on dQ. (11) 

where u x ., i — 1,2, denote the partial derivative of u w.r.t. coordinate in IR 2 . 
Without loss of generality, we assume Q is a square domain, and the mesh is uniform, 
that is, the plane IR 2 is divided into squares {(xi,^); i\h < x\ < [i\ + l)h,i 2 h < x 2 < 
(z2 + 1)^}) ii, «2 = 0, ±1, ±2, • • •, and each square is further divided into two triangles 
by a straight line x 2 = x± + ih, i integer. Take the node Xij and Qij shown as figure 
1. The elements are divided into two categories: the first category is shown as figure 2 
and the second is shown as figure 3. 

The discrete scheme for this equation form finite element method is 

u iJ+1 + u i+ ij + Uij-i + Ui-i,j - 4:Uij 

= / f(^ f u k,l ( Pk,l)'Pijdx, (12) 

where the shape function ip^ is linear continuous function defined as 

The right side of ( T2 ) is 

Isiij '■= / f(52 u k,Wk,i)Vi,jdx 

= (/ + / + / + / + / + / )f(E, u k,Ww)<Pijdx, (13) 

J A JB JC JD J E J F 
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where 



/o. 



In, , ; 



IciijE 



= / f(v>i-ij-i<Pi-ij-i + ^j-i^ij-i + Uij<Pij)<Pijdx. 

J A 

'■= / f(J2 U k,l ( Pk,l) i Pi,jdx 

= / f(ui-i,j-i(pi-i,j-i + Ui-^jcpi-^j + Uijcpi^cpijdx 
Jb 

f(^ u k,lVk,l)<Pi,jdx 



= J c f(ui-i,j(Pi-i,j + Uijipij + Uij +1 <pij +1 )<pijdx 



fQ2 u k,l<Pk,l)<Pijdx 

J f(u>i+i,j+m+i,j+i + u i,m,j + u>i,j+m,j+i)<Pi,jdx 

/ fC52 u k,Wk,i)Vi,jdx 

/ f(u i+1J+1 ip i+1 j +1 + Wij^ij + ik+ij<Pi+ij)<Pijdx 
Je 

f(J2 U k,l l fk,l)^i,jdx 

= / /(u*+ij¥>i+ij + Uijifij + Uij-iipij-Jwjdx. 

J F 



where the sub-index A, B, C, D, E, F indicate the all elements neighboring Xij in Figure 
1. 

Similar to the one-dimensional case, introducing the DEL 1-forms 

Eixiij ■■= {uij+i + u i+1J + utj-! + Ui-ij - Au i:j - I Qi . }du i:j . (14) 
The DEL condition now reads 

dE DQitj = 0, (15) 

i.e. 

du it j + i A du it j + du i+1 j A c?Ujj + duij-i A cfa^j + du^ij A dujj 

= ^{ij} + S B {i,j} + S C {ij} + S D{iJ} + S E {ij} + S F {i,j} (16) 

where, 

Sa{ij} = (J f'(u i+ljj+1 (pi +ljj+1 + Uijifij + Uij+Wij+Jpi+^+Wijdx + 

/ f'(u i+1:j+1 ip i+1:j+1 + Uijtpij + u i+1:j ip i+1:j )ip i+1:j+1 (p i:j dx) ■ du i+1J+1 A duij, 
S B{i,j} = ( / f'(u i+1:j+1 ip i+1:j+1 + Uijifij + u i:j+1 ip i:j+1 )ip i:j+1 (p itj dx + 

J ID 

/ f (ui-xjtpi-xj + Uijifij + Uij + i<pij + i)<pij +1 <pijdx) ■ du i:j+1 A du i:j , 
s c{i,j} = ( / /'(ui+ij+i^i+ij+i + Uijifij + u i+1 j(pi +1 j)(p i+1 jPijdx + 

I f'(u i+1 j(p i+1 j + uijifij + Uij-xifij-^ifi+ijifijdx) ■ du i+1 j A du^, 

J F 

So{ij} = ( / f (ui+ijifi+ij + Uij(fij + Uij-iifij-^ifij-iifijdx + 

J F 

\ f'(ui-ij-i¥i-ij-i + Uij-ifij-i + Ui,jVi,j) ( Pi,j-i<Pi,jdx) ■ duij-! A du^, 

J A 

SE{ij} = (J /'(wi-ij-i^i-ij-i + Uij-iifij-i + Uijifij^i-xj-iifijdx + 

/ /'(^i-ij-i^-ij-i + u^jifi^j + Uijif^ifi^j^ifijdx) • gK_ij_i A duij, 

J B 

S F{ij} = ( / f'iUi-lj-lfi-lj-l + Ui-ljPi-lj + Uijifij)(fi_ij(fijdx + 

J f'(ui-ijtPi-ij + Uijifij + Uij + iifij + i)fi^ijifijdx) ■ dm-xj A du^. 
Let us introduce two shift operators as 

E x f(uij) = f(u i+1 j), 
E 2 f(uij) = f(uij +1 ). 

Then the following relations can be found 

£>A{i,j} = -EiE 2 S E{i j}, 

Ss{ij} = -E 2 S D {ij } , 

s c{ij} = -E 1 S F{i j}, 

du iJ+1 Aduij = -E 2 (du i j_ 1 AdUij), 

du i+ ijAduij = -Ei(dui-ij Aduij). 



Then it is easy to check that ([16]) becomes 



D\U) Di j + D 2 T Di) j — 0. (17) 

where 

^Dij = diii-ij A duij — E 2 SE{i,j} — SF{ij}, (18) 

i~Di,j = diiij-i A du i: j — Ss{i,j} — £>E{i,j}- (19) 

It is straightforward to show that ujd and td are two symplectic 2-forms. Namely, they 
are closed w.r.t. d on the function space and non-degenerate. Then the equation (|T7|) 
is in fact the multisymplectic conservation law. Here Di and D 2 are given by 

D X = E X - 1, D 2 = E 2 - 1. (20) 

In this letter we have explored some very interesting relations between symplectic 
and multisymplectic algorithms and the simple finite element method for the bound- 
ary value problem of the semi-linear elliptic equation in one-dimensional and two- 



dimensional spaces. The details you could found in |TTJ. Although what we have found 
are certain simple boundary value problem of the semi-linear elliptic equation in lower 
dimensions and also quite simple triangulization in the finite element method, the re- 
sults still indicate that there should be very deep connections between the symplectic 
or multisymplectis algorithms and the finite element method. 

There are lots of relevant problems should to be studied and some of them are 
under investigation. 
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